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TITLE 

POWER SYSTEM 

5 DESCRIPTION 

Technical field 

The present invention relates to a power system and in particular to a method for 
voltage stabilization of an electrical power network system comprising a producing 
10 power network system side and a consuming power network side to maintain voltage. 

Background of the invention 

A power system consists of several electrical components (e.g. generators, 
transmission lines, loads) connected together, its purpose being generation, transfer 
15 and usage of electrical power. 

In a conventional On-Line Tap Changer (OLTC) the control is given by a simple 
integrator with a time delay and deadband. The size of the deadband sets the tolerance 
for long term voltage deviation. The reference signal for the integrator is the secondary 
20 voltage setpoint. This is usually kept constant at the desired secondary voltage. 

Voltage stability of a power system is defined by the IEEE Power System Engineering 
Committee as being the ability of the system to maintain voltage such that when load 
admittance is increased, load power will increase so that both power and voltage are 
25 controllable [2]. 

Voltage stability in power networks is a widely studied problem. Several voltage 
collapses resulting in system-wide black-outs made this problem of major concern in 
the power system community. 

30 

In todays state-of-the-art practice, the following methods are used to detect that the 
system is close to voltage instability: 

1. As too much power is requested by the load, the generators will start using their 
rotational energy, implying that the frequency of the voltage (50/60 Hz) will start 

35 to decrease. Detecting a low frequency has been a too slow measure to stop the 

voltage collapse in for example eastern USA in 2003. 

2. Another sign of overload is that the load voltage drops. However, it has been 
shown that neither this is a good measure for the instability of the grid. 
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Using any of the above methods (or similar), the actions taken by the power companies 
is usually one or both of the following: 

1. Connect capacitor banks, to increase the active effect that can be consumed by 
the load. If this is done in time, a voltage collapse can sometimes be avoided. A 

5 disadvantage of this method is that it makes the network more sensible to load 

variations. 

2. Disconnect certain amounts of load (load shedding). This is a very "expensive" 
measure, and therefore avoided for as long as possible by the power company. 
However this measure can prevent the whole power net from collapsing. 

10 

This invention is concerned with dynamic stability of a power systems. The inventors 
propose a dynamic feedback and feed-forward based compensation that aims at 
stabilization of the power grid. This control structure is intended to function as an 
emergency control scheme, i.e., it will be active in critical situations when the network 
15 is near voltage collapse. 

The considered power system is shown in Figure 1. It is a radial system containing a 
generator Es , a transmission line with impedance Z in , a transformer with an on-line 
tap changer (OLTC) and a load with impedance Zld . The on-line tap changer regulates 
20 the voltage on the load side at a desired value V r£f . The load itself dynamically 

changes its impedance. Most of the loads are such that they try to absorb a certain 
amount of power. That implies that when the load voltage drops, the loads will 
decrease their impedance to keep power constant. 

25 There are two control loops in this system, acting independently of each other. 

• The On-Line Tap Changer (OLTC) in the transformer, which tries to keep the 
voltage on the load side constant at the reference value V ref , 

• The load itself can be viewed as a control system, which changes its impedance 
(or equally admittance) in order to absorb a given power. 

30 

The problem is that these two independent control loops can, due to their non-linear 
interaction, drive the system to voltage instability even if the system could handle the 
power required by the load. 

35 Summary of the invention 

This work proposes a general method that momentarily changes the behavior of the 
OLTC when the line and/or load impedance changes such that the system is driven into 
the critical operation regime. 
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It is important to again point out that the proposed control structure is meant to 
operate in case of dynamic instabilities. This means that after a line and/or load 
impedance change (for example due to a line failure or an increase of power request 
5 from the load) the power grid is still statically capable of transferring the load power 
request. 

In particular the method of the invention is characterized in that the power transfer Y L d/ 
wherein Y LD is power load impedance, is dynamically maintained below the loci for 
10 maximum power transfer, n 2 Y UD Z LN = 1, wherein Y LD is power load impedance, Z LN is 
transmission line impedance and n is transformer ratio, preferably Y L d is maintained at 
a stable equilibrium. 

The present invention makes use of a mathematical model: 

15 

For ease of reference a list of used variables is compiled below: 



Zld = Zld e J ® - load impedance, 
Yld = II Zld - load admittance, 



20 



Zln = Zlnc j -transmission line impedance, 
E s = E s e j0 - generator voltage, 



Vi- voltage on the primary side of the transformer, 



V 2- voltage on the secondary side of the transformer, 
rc-transformator ratio, 
y r ef - reference voltage, 



25 



Ii~ current in the primary winding of the transformer, 



/ 2- current in the secondary winding of the transformer 



For the system in Figure 1, some basic relations can be stated [4]: 



30 



y 2 /Vi = /i//2 = n 



Es = IlZ ln+ V 1 = l2(nZ In+l/nZLD) 




cos<3> 
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The function is a nonlinear function that determines the typical dependence of the 
active power on the line and load impedance (Figure 2). Initially, for increasing Y LD r 
the active power will increase. However, after a certain load admittance the transfered 
5 active power starts to decrease. For Zld/h 2 = Z in a maximum active power will be 
transmitted through the line. 

Then for a constant active power load, a suitable model is: 

10 YLD = Pref~PR = Pr ef - Es* ■ ~ , , COS O (1) 

|Zln+ZLD/tt 2 | 

while the OLTC can be approximated by an integrator: 

_ Zld/h 

n = Vref-Es (2) 

I Z in+ Zld I n I 

In order to understand the behavior of the proposed model, consider first the dynamical 
15 system in equation (1). Due to the built-in non-linearity, the system can have two 

equilibrium points corresponding the reference active power (see Figure 2). It can be 
shown that the one to the left of the peak is stable while the other is unstable. This will 
determine the typical behavior of a power system. After achieving the maximum value 
of the transfered active power, if the load admittance continues to increase, the system 
20 enters the unstable region. This will lead to instability if the load admittance achieves 
the value corresponding to the unstable equilibrium point. 

Simulation results for the above model are shown in Figure 3. The variable in the plot 
are the maximum transferable active power, the transfered active power and load 

25 impedance. In this scenario the load is trying to absorb an active power of 0.7 (dashed 
line). The initial value for the line impedance is 1. At t-7S a fault is simulated in the 
line by changing its impedance to 1.5. As shown in the first sub-plot, this implies that 
the maximum power that can be transferred through the line will drop just below 0.7. 
The load tries to absorb the desired active power by reducing its impedance (see the 

30 second and third sub-plot). However since that power is not achievable, the system will 
end up in instability and voltage collapse. 

Considering both equations (1) and (2) in the model, similar qualitative behavior is 
retain as for the scalar case. Figure 4 shows the vector field near the equilibrium points 

35 (marked with asterisks). The dashed line is given by the curven 2 7z,DZin = 1, i.e. the loci 
of maximum power transfer (this happens if the line impedance and the load 
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impedance are equal). Notice the unstable behavior to the right of this curve. 

The present mathematical model is able to capture two instability scenarios. 

1. The first case is shown in Figure 3, where due to some fault in the transmission 
5 line the system is no longer able to transfer the requested active power. This 

corresponds to the situation when the system has no real equilibrium points. This 
is the classical case, which can be analyzed even with static methods. 

2. Another instability scenario is when a stable equilibrium point exists, but where 
the system ends up in instability due to some transients. In Figure 6, at 50 time 

10 units, a fault in the transmission line is simulated by a step increase of the line 

impedance. This step is such that a stable equilibrium point still exists, that is, 
the network should be able to transfer the requested active power. However, due 
to the fact that the operating point is close to the maximum transferable active 

power, an overshoot in Yn , will drive the system in the unstable region and the 
15 voltage will collapse. 

The methods described in this paper adds stability margins so that the risk of the 
second scenario is significantly reduced. The stabilizing property of the methods will 
also help restoring stability after an overload condition when load shedding has been 
20 applied. 

The proposed methods comes in before the methods 1 and 2 above would be applied. 
This way, adds no inconvenience to the customers while preserving stability. If stability 
cannot be maintained in spite of these methods (due to too large power demands), the 
25 methods above should be applied. 

As can be seen in Figure 4, it is desirable to move the system away from the unstable 
region above the stability limit (dashed curve). Since the load dynamics cannot be 
changed (except by load shedding), we suggest to momentarily alter the transformer 
30 ratio n so as to avoid the unstable region. 

The following sections describe how this can be done in practice, indirectly, by changing 
the voltage reference V ref given to the standard OLTC. 

35 A block diagram over the structure of the proposed compensator is shown in Figure 7. 

The compensator consists of two susbsystems. The first susbsystems consists of a 
feed-forward compensator and the second consists of a feedback controller. 
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The goal of the feed-forward compensation is to improve the convergence ratio of the 
system in case of a fault in the transmission line. In other words, the compensator will 
drive the system to the stable equilibrium point in case of a line fault. However, this 
method works only if, after the fault the system is still the stable region (i.e. 

5 n 2 YLDZht<l). 

The idea of using such compensation is suggested by the structure of the presented 
simplified model. It is rather straightforward to show that the line impedance Z In ac t s 
as a load disturbance on the system, similarly to ^ r #. In addition, the line impedance 

10 can be considered measurable. It is natural then to use a feed-forward compensation 
from the line impedance in order to diminish the influence of line faults. If the 
transformer ratio n would be directly accessible for control purposes, the transient 
influence of line fault could be (at least theoretically) completely removed. Although 
only V r ef is accessible, it is still possible to considerably improve the line-fault behavior 

15 of the system. 

This compensating subsystem aims to prevent the grid from entering an unstable 
operating regime. For this it uses information about the line impedance. 

20 A suitable feedforward compensation is given by the first order filter 



where T, T d are tuning parameters. 

25 In case the system enters the unstable region (i.e. h 2 YldZ in > 1 ), another control 
strategy has to be applied, which is described in the next section. 

When the system is in the unstable region, it is desirable to drive it back to the stable 
30 operation regime. This can be done by reducing the reference voltage as long as the 
system is in the unstable region, Such a compensation can be achieved by a static 
nonlinear feedback. In Figure 4, as a result of the compensation, the vector field above 
the line y^YldZ^ = lwill point inwards (see Figure 5). It can be seen in the the plots 
that the region of attraction for the stable equilibrium point has been considerably 
35 increased. 



It is to be mentioned here that the idea of using the distance from the peak of the 
function fin voltage stability studies has been recently proposed in [3]. However, it has 
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never been used (to the best of the authors knowledge) for dynamic compensation of 
the voltage reference signal. 

Thus the second control subsystem aims to drive the grid from the unstable operation 
regime to the stable operation regime. For this it uses information about the line 
impedance, load impedance, and transformer ratio. 

A suitable feedback controller is: 

^ = -max(0,a(n 2 F LD -1/Z h )) 

where a is a tuning parameter that is influencing the region of atraction of the 
equilibrium point. 



In order to obtain more realistic simulation results the initial design model has been 
modified as follows: 

• the dynamics have been scaled according to the benchmark model [5], 

• additional dynamics have been introduced for the load argument, cp, 

• load shedding input k has been added, 

• saturation and quantization is introduced on the transformer ration n. The latter 
is intended to simulate the mechanical tap-changer, 

• since the tap-changer is inherently a discrete system, a discrete time 
representation of the OLTC dynamics is used. Notice that the tap-changer can 
make only one step at the time. 

• in order to avoid chattering, an OLTC system usually contains a dead-zone on the 
control error. 

This way the simulation model is the following: 



Y = l/T 



\l - k) Pr ef - Es 2 , ~ Zl £ /n2 2 - 2 COS ^ 



<S> = {l~k)Qref-\IT<&-E s 2 ^ *% , 2)2 sinO 



Zld/u 2 
\Zm+ZLD/n 2 | 

7]{t + h) = Tj(t) + qsign(e(t)) 

e(t) = dzn\ Vref - Es ■ „ Zl Z ,h - 0 ■ 
^ \Zm+ZLD/n 2 | j 

n = sat{rf) 



The saturation on n has the limits 0.75 , H max — 1.25 , and the dead-zone has the 

limits ±0.03. The chosen quantization step q is 0.027. The chosen sampling time is 30 
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seconds, which approximates the mechanical delay of the tap-changer and the OLTC 
delay timer. 

The three-stage control system consists of the following compensator: 

30s 

5 • feed-forward compensation: Hff(s) = has a "dirty-derivative" character 

20S-M 

with the low-pass filter having its time constant comparable with that of the 
controlled system. 
• feedback compensation: 

Vjb = - max(0, a(n 2 YLD -II Z in)) . The parameter a influences the region of attraction 
10 of the equilibrium point. In the simulations a=l.l. 

The first two control signals ( and ) augment the reference value as follows: 

e(t) = dzn(Vre f + V ff + Vjt - & , - Z *£ / " 2 , ) 

15 

where dzn is the dead-zone function. 

However, a more complex augmentation is also possible, e.g. V ff \s conditioned by 

20 In the simulations, the following parameters have been used: 

Vre/ = 1.1, Pre/ = 0.78, Es = 1.5,7"= 60, and 6=1.47 radians. In addition, in the first 
simulation scenario (Figure 8) the reference reactive power \sQref = 0.16. The scenario 

consists of a line tripping at £=800 seconds, when the line impedance Z ln j s increased 
from 1 to 1.2. The first 800 seconds in the simulations represent the initial transient to 

25 the studied equilibrium point and it has no physical interpretation. At the moment of 

the fault, Vff shows a significant increase. However, since the new equilibrium point is 
not achieved the system ends up in the unstable operating region (at around 1100 
seconds). This will trigger the second stage of the controller, decreasing V ^ . This will 
result in a decrease of the overall voltage reference value such that the system is 

30 brought back in the stable region. Notice that throughout the entire control sequence, 
the third control stage (load shedding) is not engaged, i.e. k=0. 

It is important to remark that the first step (i.e. V ff ) j S sensitive to the fault timing 
due to the low sampling frequency. Similarly if multiple steps (e.g. two) would be 
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possible, the performance would increase significantly. Nevertheless, even in the case 
of the state-of-the-art OLTCs, where the delay timer is inverse proportional to the 
control error, considerable improvements can be obtained in compensating for line 
tripping. 
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